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En este artículo se presenta un nuevo procedimiento para la 
regionalización de los escurrimientos máximos instantáneos 
registrados en estaciones hidrométricas de las 37 regiones 
hidrológicas en la que se encuentra dividida la república 
mexicana. Las estaciones hidrométricas seleccionadas 
se caracterizan por medir escurrimientos de cuencas no 
influenciadas por la presencia de embalses o lagos de 
regulación. La relativa simplicidad del método permite 
aplicarlo en el caso de cuencas con escasa información 
hidrométrica, pero que cuenten con datos de precipitaciones 
máximas anuales, y con características fisiográficas y 
cinemáticas medibles, entre ellas área de la cuenca y tiempo 
de concentración, además de retención potencial máxima 
para caracterizar pérdidas. Se utiliza la media de los gastos 
máximos anuales (MOMIA) como concepto eje, de tal manera 
que al relacionarlo con características medibles en cualquier 
cuenca (área, tiempo de concentración, promedio espacial 
de la media de precipitaciones máximas anuales, retención 
potencial máxima), es posible estimarla de forma sencilla, 
aun para cuencas no aforadas; por otro lado, se agrupan 
estaciones homogéneas en cuanto al comportamiento 
estadístico de los gastos máximos anuales, para estimar 
factores que permiten pasar del valor de la MOMIA al 
correspondiente a distintos periodos de retorno. 

Palabras clave: regionalización, algoritmos genéticos, 
escurrimiento superficial, media de los gastos máximos 
anuales. 
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G. E. (September-October, 2016). Regional Models for Instant 
Maximum Annual Flows in Mexican Republic. Water Technology 
and Sciences (in Spanish), 7(5), 15-32. 


A new procedure for maximum annual floods regionalization in 
hydrometrical stations belonging to 37 hydrologic regions of the 
Mexican Republic is introduced in this study. The selected stations 
measures runoff of basins free of reservoirs or regulation lakes 
upstream. The relative method's simplicity allows its application 
in case of ungauged basins but they must have rainfall records and 
measurable physiographic and kinematic characteristics, for example 
the drained area, concentration time, besides the maximum potential 
soil moisture to consider the lost. The annual maximum mean inflows 
(MQMIA) were used as axis concept, such that relating it with 
measured characteristics in any basin (drained area, concentration 
time, annual maximum spatial mean rainfall, maximum potential 
soil moisture) to estimate MOMIA in a simple way is possible, even 
for ungauged basins; on the other hand homogeneous stations are 
clustered taking into account the statistical behavior of the annual 
maximum runoff, to get factors that transforms the value of the 
historical MOMIA to that one corresponding to different return 
periods . 


Keywords: Regionalization, genetic algorithms, surface runoff, 
annual maximum mean inflows. 
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Introducción 


En la república mexicana, la información cli- 
matológica es considerablemente mayor que 
la hidrométrica. En 2012, México contaba con 
3 817 estaciones climatológicas oficiales ope- 
radas por la Comisión Nacional del Agua (Co- 
nagua) y la Comisión Federal de Electricidad 
(CFE), mientras que las estaciones hidrométricas 
con registro son 2 226 de las cuales se selec- 
cionaron 309 con 20 o más años de registro, y 
cuyo funcionamiento no está alterado con obras 
hidráulicas (Conagua, 2014a). 

La estimación de valores máximos del escu- 
rrimiento con fines de estudios hidrológicos es 
un importante problema por resolver. El análisis 
regional permite hacer agrupaciones de sitios 
con comportamiento hidrológico estadística- 
mente similar (escurrimientos, precipitaciones), 
con lo que se logra la construcción de registros 
lo más grandes posibles para obtener una fun- 
ción de distribución de probabilidades de toda 
la zona considerada (Domínguez et al., 2012). 
Por otra parte, al relacionar estos resultados con 
las características fisiográficas y de precipitación 
de las cuencas, es posible deducir datos en sitios 
con escasa o nula información a partir de los 
obtenidos en los sitios aforados (Domínguez, 
Arganis, Gonzáles, Carrizosa, £ Guzmán, 2013). 

En distintas partes del mundo se han llevado 
a cabo estudios de regionalización basados en 
diversas técnicas para agrupar cuencas, toman- 
do en cuenta datos medidos, así como caracte- 
rísticas fisiográficas, y datos de uso de suelo y 
climatológicos; entre ellos resaltan los trabajos 
de Domínguez (1983), quien propuso modifica- 
ciones a las ecuaciones de estimación del gasto 
a partir de características fisiográficas del Flood 
Studies Report (1975), documento valorado tam- 
bién por Hosking, Wallis y Wood (1985); por otra 
parte, Gottschalk (1985), y Nathan y McMahon 
(1990), identificaron cuencas regionalmente ho- 
mogéneas en Suecia y Australia, respectivamen- 
te, a partir del análisis de los escurrimientos, 
usando estudios de componentes principales y 
regresiones no lineales múltiples; otros modelos 
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de regionalización basados en técnicas estadísti- 
cas y relaciones de tipo no lineales potenciales, 
análisis multivariado y manejos de sistemas de 
información geográfica, así como validación 
de metodologías de regionalización son los de 
Domínguez (1981); Carrizosa (1997); Allasia y 
Avruch (2000); Erazo (2004); Cabrera, Guardado, 
Peláez y Gonzáles (2004); Azevedo, Demetrius 
y Reis (2008); Zamanillo (2008); Isik y Singh 
(2008); Ganancias (2010); Olguín, Rivera y Sa- 
las (2011); Luna y Domínguez (2013), y Byzedi, 
Saghafian, Mohammandi y Siosemarde (2014). 

En este trabajo se presenta un procedimiento 
para obtener funciones que determinan la media 
del gasto máximo instantáneo anual (MOMIA) 
para las cuencas de una región, en función 
de sus características fisiográficas fácilmente 
medibles y de datos de precipitación. Estos re- 
sultados pueden ser aplicables después en sitios 
de la región que no cuenten con información de 
escurrimientos, pero que sí se pueden describir 
a partir de su fisiografía y precipitación. 

Además, tomando en cuenta la similitud 
en las características estadísticas de los gastos 
máximos registrados en grupos de cuencas, 
se desarrolló un análisis regional que permite 
determinar eventos asociados con distintos pe- 
riodos de retorno mediante factores regionales 
que multiplican a la MOMIA de las estaciones 
hidrométricas de cada grupo. 


Metodología 


Extracción de gastos instantáneos y de gastos 
medios diarios 


La fuente oficial de datos medidos por las esta- 
ciones hidrométricas en México se encuentra en 
el Banco Nacional de Datos de Aguas Superfi- 
ciales (BANDAS) de la Comisión Nacional del 
Agua y del Instituto Mexicano de Tecnología 
del Agua (IMTA) (Conagua, 2014a); dicha base 
de datos cuenta con archivos con nombre “con- 
cepto regiónclavedeestación.dbf” de los que se 
extraen los registros para los gastos máximos 
instantáneos (concepto HD) y para los gastos 
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medios diarios (concepto DD). Se elaboró un 
programa codificado en Visual BasicO para la 
conversión de estos archivos del BANDAS a 
archivos de un hoja de cálculo ExcelO, y que 
además realiza el cálculo de los gastos medios 
máximos para distintas duraciones y los esta- 
dísticos media, desviación estándar, variancia y 
coeficiente de asimetría. A dichos archivos se les 
inserta el registro de gastos máximos instantá- 
neos, también obtenido con un programa desde 
los archivos del BANDAS. Después se realizan 
comparaciones entre los valores de los gastos 
máximos instantáneos con respecto a los gastos 
medios máximos de un día, y en los casos que 
se Observen mayores los medios diarios que 
los instantáneos, se opta por eliminar el dato o 
bien se revisa la información de gastos máximos 
mensuales (concepto DM), que en muchas oca- 
siones incluye valores máximos no incluidos en 
el concepto HD. Por último, cuando subsiste la 
duda y el valor correspondiente puede ser im- 
portante, se recurre a conseguir la información 
original. 


Cálculo de estadísticos de los gastos máximos 
instantáneos 


Se determinan la media, variancia, desviación 
estándar y coeficiente de variación, además del 
valor máximo y mínimo del registro de gastos 
máximos instantáneos. 


Cálculo de los atributos correspondientes a 
cada cuenca 


Area drenada 


El valor del área drenada se obtiene de la 
descripción de estaciones hidrométricas del 
BANDAS de la Conagua por ser la fuente oficial 
de la nación. En los casos en que no se cuenta 
con la información, el área drenada (figura 1) 
se calcula a partir de una herramienta digital 
del tipo ArcmapO o ArcviewO, ubicando a la 
estación hidrométrica con sus coordenadas geo- 
gráficas y utilizando la extensión GeoHMSO. En 
este estudio se calcularon las áreas drenadas de 


Figura 1. Área drenada. Ejemplo de la estación hidrométrica 18434 El Fraile. 
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todas las estaciones consideradas y el resultado 
se verificó con el que proporciona BANDAS, de 
manera que si se obtenían diferencias impor- 
tantes, se utilizaron los valores obtenidos con 
ArcviewO. 


Tiempo de concentración 


El tiempo de concentración (tc) se determina con 
la longitud (L) y pendiente del cauce principal 
(S), con la ecuación experimental de Kirpich 
(Aparicio, 2012): 


po” 
tc =0.000325=— (1) 


90.385 
Donde tc en h, L en m y S en decimales es 
determinada con el método de la pendiente 
media. 


Precipitación media hp 


Se identifican las estaciones climatológicas ubi- 
cadas dentro o cerca de la cuenca asociada con 
cada estación hidrométrica analizada; se calcula 
la media de las precipitaciones diarias máximas 
anuales para cada estación climatológica y se 
obtiene su promedio espacial. 


Número de escurrimiento N y retención potencial 
máxima S 


El número N se determina con el tipo y uso 
de suelo a partir de los cuadros del US Soil 
Conservation Service. A partir del número de 
escurrimiento se obtiene el espesor de suelo S 
(infiltración potencial máxima). En Domínguez 
et al. (2008) se puede encontrar en forma detalla- 
da la determinación de S, dada por la ecuación 


(2): 


g 2540 =((25.4x N) 
N 


(2) 


Donde N es el número de curva y S es el 
espesor de suelo que corresponde a la retención 
potencial máxima, en cm. 

De la ecuación (2) se observa que el número 
de curva N depende de un solo parámetro y 
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de acuerdo con Paz et al. (2010), esta particu- 
laridad lo hace atractivo para usar N como un 
método hidrológico operacional; en este estudio 
se utilizó la retención potencial máxima como 
parámetro de caracterización de las pérdidas al 
escurrimiento. 


Agrupación de regiones 


Se consideran primero cada una de las 37 re- 
giones hidrológicas señaladas por la Conagua 
(figura 2). 

Enseguida, se construyen gráficas de las 
áreas drenadas contra la media de los gastos 
máximos instantáneos anuales MOMIA de las 
estaciones hidrométricas de cada región hidro- 
lógica (o grupo de regiones hidrológicas) para 
establecer grupos de estaciones, con base en si 
sus valores de MOMIA están por arriba, por 
abajo o casi sobre la curva de ajuste (figura 3). 

También es de utilidad para verificar la 
agrupación, analizar la ubicación espacial de las 
estaciones en esquemas topológicos, como el de 
la figura 4, a los que se les añade el coeficiente 
de variación, para observar la homogeneidad de 
los grupos. En el caso de la región 20, la figura 
muestra que la agrupación obedece al tipo de 
escurrimiento: si es del grupo 1, se trata de una 
corriente principal; el grupo 2 corresponde a 
ramales secundarios, y el grupo 3 a ramales 
terciarios. 


Registros normalizados para cada grupo de 
estaciones hidrométricas 


Los datos anuales del gasto máximo de cada 
estación hidrométrica se dividen entre su media 
y por cada grupo de estaciones hidrométricas 
se forma un solo registro con todos los datos 
normalizados, el cual se trabaja como una serie 
anual. 


Análisis estadístico de series normalizadas 


La serie anual formada por los datos normaliza- 
dos de cada grupo de hidrométricas se analiza 
estadísticamente, ordenándolos de mayor a 
menor; se obtiene su periodo de retorno con 
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Figura 2. División del país en 37 regiones hidrológicas. Fuente: Conagua, 2014b. 
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Figura 3. Ejemplo de gráfica de la MOMIA us. el área drenada (RH18). 


una ecuación empírica (en este caso se usó la mejor ajuste. Con la función de distribución de 
de Weibull; Aparicio, 2012) y se hace el ajuste mejor ajuste se calculan valores de los gastos 
con distintas funciones de distribución, de las máximos instantáneos normalizados para dis- 
cuales se selecciona aquella que proporciona el tintos periodos de retorno. 
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Figura 4. Ubicación de los grupos de hidrométricas. Ejemplo de una porción de la RH20. 


Los valores así obtenidos para cada grupo 
de estaciones hidrométricas constituyen factores 
normalizados, de tal manera que al multiplicar- 
los por la MOMIA de cada estación hidromé- 
trica del grupo correspondiente se obtiene una 
estimación regional de gastos máximos para 
distintos periodos de retorno. 

Debido a que las muestras de cada grupo 
contienen muchos más elementos que cuando 
se trabaja con cada estación por separado, los 
resultados del análisis descrito permiten una 
estimación confiable de los gastos de diseño 
para las estaciones hidrométricas usadas en 
el análisis o para algunas otras que no fueron 
utilizadas porque su registro es corto, pero que 
puede considerarse suficiente para estimar la 
media (p. ej., cuando dicho registro es de entre 
15 y 20 años). En los sitios sin registro histórico 
se requiere un procedimiento alternativo para 
estimar la media, como el que se describe a 
continuación. 


Curvas regionales para estimar la MOMIA 


Se proponen modelos no lineales de tipo 
potencial como los de las ecuaciones (3), (4) 
y (5): 
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MOMIA =C14? (3) 
MOMIA =C14 tc (4) 
MOMIA =C1V “tc st (5) 


Donde MOMIA es la media de los gastos 
máximos anuales obtenida para cada estación 
hidrométrica de la región o grupo de regiones 
considerada, en m?*/s; A, el área de la cuenca, 
en km?,; tc, el tiempo de concentración en h; 
V, un indicador del volumen de precipitación 
en la cuenca, obtenido al multiplicar el área 
(km?) por la precipitación hp (mm) (determi- 
nada como el promedio en la cuenca de la 
media de las precipitaciones diarias máximas 
anuales. El valor correspondiente a dicha 
precipitación, así como los valores asociados 
con distintos periodos de retorno y duración 
pueden consultarse en Domínguez, Carrizo- 
sa, Fuentes, €: Galván, 2010, en la que se usó 
la información de poco más de dos mil esta- 
ciones pluviométricas con más de 20 años de 
registro), y S es la retención potencial máxima 
del suelo en cm. 

Se determinan los parámetros C1, (2, C3 
y C4 de los modelos (3), (4) y (5); para ello se 
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utilizaron algoritmos genéticos que son herra- 
mientas del cómputo evolutivo (Holland, 1975; 
Goldberg, 1989; Jiménez, 2004; Domínguez, 
Fuentes, Arganis, £ Mendoza, 2009; Rodríguez, 
Arganis, Cruickshank, £ Domínguez, 2012). 


Selección del modelo 


La selección del modelo debe hacerse tomando 
en cuenta la información disponible en la cuenca 
analizada, los coeficientes de correlación entre 
los datos medidos y calculados por el modelo, y 


los coeficientes de correlación se construyen grá- 
ficas de los datos medidos contra los calculados 
por cada modelo en relación con una función 
identidad (figura 5). 

En la figura 5, el modelo que presenta un ma- 
yor coeficiente de correlación y determinación 
(r =0.9131 y R?= 0.8338) es el que relaciona a la 
MOMIA con el volumen, tiempo de concentra- 
ción y espesor del suelo; este modelo considera 
de forma implícita al área de la cuenca y preci- 
pitación, además de procesos de infiltración y 
tiempo de respuesta de la cuenca al volumen 
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Figura 5. Ejemplo de resultados de ecuaciones regionales de la media calculada en relación con la media medida y la función 
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Aplicación y resultados 


Estimación de factores regionales para 
cuencas aforadas 


Se escogieron 309 estaciones hidrométricas 
con más de 20 años de registro histórico, con 
la revisión ya descrita. Se procuró que dichas 
estaciones no estuvieran influenciadas por 
embalses aguas arriba. Los grupos de regiones 
hidrológicas se establecieron como se indica en 
el cuadro 1 (ver figura 2). 

Para dar como ejemplo el caso de la región 
hidrológica 18 (RH18), con el procedimiento 
descrito se comparó la MOMIA contra el área 
de la cuenca (figura 3) y con base en esta figura 
se establecieron tres grupos de hidrométricas 
(columna 9 del cuadro 2). El grupo 1 corres- 
pondió a las estaciones que quedaban por arriba 
de la curva; el grupo 2, a las hidrométricas que 
quedaron por debajo de esa curva; el grupo 3, 
a los datos que quedaron aproximadamente en 
la curva o muy cercanas a ella. 


Cuadro 1. Agrupación de regiones. 


Grupos RH 
1 1y3 
2 7y8 
3 9 
4 10 
5 11 
6 12 
Z, l13a16y19 
8 18 
9 20 y 21 
10 220, 
1 23 
12 24 
13 25 
14 26 
15 27 
16 28 y 29 
17 30 
18 36 y 37 
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Además, en el cuadro 2 se indican los 
estadísticos del gasto máximo instantáneo, en 
particular la media y el área drenada por cada 
hidrométrica. 

Por otro lado, se normalizó la información 
de los gastos máximos anuales de las estaciones 
correspondientes a cada grupo y se construyó 
un gran registro por grupo, al que se le hizo un 
análisis estadístico para determinar los gastos 
máximos instantáneos normalizados para 
distintos periodos de retorno (en la figura 6 se 
muestra la extrapolación para el grupo 1 (G1) 
y en el cuadro 3 se presenta el resumen de los 
factores extrapolados para los tres grupos (Gl, 
G2 y G3) de la RH18). 

En el cuadro 4 se muestran los factores cal- 
culados para todos los grupos. 

Los valores del cuadro 4 se utilizan como 
factores, de manera que al multiplicarlos por la 
MOMIA de cada estación se obtienen las esti- 
maciones de sus gastos máximos asociados con 
distintos periodos de retorno, como se muestra 
en el siguiente ejemplo. 


Ejemplo del uso de factores para cuencas 
aforadas 


El siguiente ejemplo tiene la finalidad de mos- 
trar la aplicación de los resultados del cuadro 
4, así como contrastar los resultados que se 
obtienen del análisis regional con los que se 
obtendrían analizando de forma individual los 
gastos máximos anuales de una cuenca aforada. 

Si se sabe que la MOMIA histórica de la 
estación 18481 La Caimanera es de 2 643.804 
m'/s, y que pertenece al grupo 3 de la RH18, 
considerando el factor de 3.09 del cuadro 4, el 
gasto máximo instantáneo de diseño —corres- 
pondiente a un periodo de retorno de 50 años, 
considerando un análisis regional, en el que se 
toman en cuenta los datos de ocho estaciones— 
es de 2 643.804 x 3.09 = 8 166.99 m?/s. Si se hace 
el análisis estadístico de los gastos máximos 
registrados en esa estación, sin regionalización, 
el gasto de diseño resulta de 11 008.09 m*/s, 
obtenido con una función doble Gumbel. Es 
decir, la regionalización da un resultado menos 
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Cuadro 2. Grupos RH 18. Estadísticos del gasto máximo instantáneo, área drenada y agrupación de hidrométricas. 


Estación Nombre Máximo Media Desv. est. | Coef. var. | Mínimo Área, km? | Grupos 
18481 La Caimanera 9906 2 643.804 | 1737.300 0.657 826 80 592.475 3 
18230 Santo Tomás 3881 1 569.484 748.473 0.477 405.21 59 153 2 
18439 San Juan Tetelcingo 3264 1 329.194 572.351 0.431 382 45 361.7 2 
18434 El Fraile 1569.74 498.872 274.995 0.551 201 ias 2 
18494 Los Panches 4778 1167.623 | 1009.516 0.865 128 11 596.4 1 
18459 Río Grande 136 73.696 23.224 0.315 10 943.5 3 
18432 Ixcamilpa 3210 586.302 472.631 0.806 159 4 946.6 3 
18487 Los Pinzanes 5 422 924.308 896.867 0.970 146 5 525.1 1 
18195 Ziritzícuaro 1 230.2 533.451 296.523 2.440 155 1 701.4 1 
18201 El Cajón 1116 264.983 242.082 0.914 28.12 395 3 
18495 La Pastoría 2873 891.102 527.682 0.592 365 2 655.9 1 
18461 Zitácuaro 297.47 112.922 56.465 0.500 40.9 372.5 3 
18129 Los Granjenos 258 50.957 50.421 0.989 11 171.4 3 
18281 La Limonera 673.56 316.476 114.098 0.361 173 1 681.3 3 
18460 Río Chiquito 875 71.904 115.542 1.607 12 247.8 3 

Doble Gumbel. P = 0.92 Tr = años 
e 
EE O 
10 
9 
8 - 
7 
2 6 $ 
BO5 . 
a, = _e a A _ _ 
3 ¿at : - 
2! ' 4 - A 
1 
0+ j Y r ' 
2 1 0 1 2 3 4 5 6 Z 8 9 10 
Variable reducida, Z =-In [In (Tr /(Tr-1) ) ] + Medidos e Calculados 


Figura 6. Ejemplo de la extrapolación probabilística para el Grupo 1 (G1) de la Región Hidrológica 18. 


conservador que tomar en cuenta la cuenca en 
forma independiente de los sitios restantes. 

Cabe mencionar que los factores del cuadro 
4 pueden aplicarse para cuencas “vírgenes”, 
aunque cuenten con un número de años de 
registro menor (p. ej., 15 años). 


Estimación de la MQMIA en sitios no 
aforados 


Por otra parte, con el área de la cuenca y datos 
de precipitación, infiltración y cinemáticos 
(tiempo de concentración), se obtuvieron mo- 
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Cuadro 3. Factores de diseño normalizados. 


e RH 18 
Gl G2 G3 
Años 
0.82 0.89 0.86 
1.34 1.31 1.27 
10 1.83 1.67 1.63 
20 2.56 2.02 2.14 
50 3.78 2.45 3.09 
100 4.64 2.75 3.82 
200 5.47 3.05 4.54 
500 6.54 3.43 5.45 
1000 7.34 3.73 6.14 
2 000 EI 4.02 6.83 
5 000 9.18 4.4 7.74 
10 000 9.96 4.7 8.41 
Cuadro 4. Factores de diseño normalizado para todos los grupos de regiones. 
= RH1a3 RH 1 (sin 3) RH7y8 RH9 RH 10 
G1 G2 G3 Gl G2 G3 G1 Gl G2 G1 G2 
Años 
2 0.32 0.49 0.53 0.32 0.42 0.33 0.73 0.83 0.77 0.68 0.88 
5 1.03 1.14 1.34 1.03 0.92 0.94 1.34 1.34 1.37 1.32 151 
10 4.16 3 2.32 4.16 3.89 3.85 1.91 1.87 2.03 2.22 1.92 
20 5.89 5.57 5.19 5.89 6.12 7.86 3.18 2.45 2.87 3.42 2.32 
50 7.64 8.34 8.85 7.64 8.37 11.91 6.07 3.04 3.84 4.72 2.83 
100 8.86 10.29 11.29 8.86 9.93 14.73 7.99 3.43 4.51 5.61 3.22 
200 10.04 12.2 13.62 10.04 11.45 17.47 9.82 3.82 5.15 6.48 3.6 
500 11.58 14.68 16.64 11.58 13.43 21.04 12.16 4.31 5.99 7.6 4.1 
A 1000 12.74 16.55 18.91 12.74 14.93 23.71 13.92 4.68 6.62 8.44 4.49 
2 2 000 13.9 18.43 21.13 13.9 16.42 26.43 15.66 5.05 EZ 9.27 4.87 
E 5 000 15.4 20.84 24.11 15.4 18.34 29.97 17.91 5.55 8.06 10.39 5.37 
El 10 000 16.54 22.68 26.33 16.54 19.89 32.61 19.62 592 8.68 11.27 5.75 
E 
E D RH 11 RH 12 RH 13 a 16 y19 RH 18 
E G1 G2 G3 G1 G2 G3 G1 G2 G3 G1 G2 G3 
E 2 0.7 0.68 0.8 0.87 0.75 0.8 0.69 0.72 0.77 0.82 0.89 0.86 
ps 5 1.37 1.67 15 1.38 1.37 1.39 1.41 1.3 14 1.34 1.31 1.27 
dE 10 2.27 2.48 2.06 1.83 2.03 2.01 2.49 1.83 1.97 1.83 1.67 1.63 
E 20 3.12 3.17 2.77 2.34 Sus 2.79 32 2.88 2.81 2.56 2.02 2.14 
S 50 4.12 4.02 3.97 3 4.43 3.67 3.98 5.55 4.13 3.78 2.45 3.09 
7 100 4.84 4.63 4.95 3.48 5.26 4.26 4.52 7.39 5.08 4.64 E 3.82 
e 200 5.54 5.23 5.93 3.95 6.05 4.82 5.06 9.14 6 5.47 3.05 4.54 
E 500 6.46 6.02 pz 4.56 7.07 5.55 5.76 11.4 7.18 6.54 3.43 5.45 
$ 1000 7.16 6.62 8.15 5.02 7.83 6.1 6.28 13.09 8.06 7.34 3.73 6.14 
E 2 000 7.86 7.21 9.11 5.47 8.58 6.65 6.81 14.77 8.94 8.12 4.02 6.83 
S 5 000 8.75 8.01 10.34 6.07 9.58 7.36 75 16.93 10.1 9.18 4.4 7.74 
S 10000 | 9.48 8.6 11.33 6.55 10.38 7.89 8.02 18.68 10.97 9.96 4.7 8.41 
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Cuadro 4 (continuación). Factores de diseño normalizado para todos los grupos de regiones. 


RH 20 y 21 RH 22 RH23 RH24 RH 25 
id G1 G2 G3 G1 G2 G1 G1 G2 G3 G1 G2 
Años 
2 0.79 0.82 0.86 0.65 0.7 0.84 0.65 0.54 0.74 0.81 0.57 
5 1.29 1.43 1.34 1.43 1.42 1.41 11.577 1182 1.34 1.42 1.35 
10 1.93 2.06 1.77 2.17 2.18 1.89 2.31 2.05 1.85 2.02 1.99 
20 2.92 2.66 2.26 3.44 3.45 2.41 3.73 3.81 2.74 2.76 2.43 
50 3.57 3.25 2.9 5.53 4.95 3.09 5:22 8.45 5.45 3.69 2.93 
100 3.99 3.65 3.37 6.96 5.94 3.59 6.22 11.51 Aoi 4.33 3.29 
200 4.38 4.04 3.83 8.33 6.87 4.08 7.18 14.39 9.47 4.95 3.65 
500 4.89 4.54 4.43 10.09 8.08 4.73 8.43 18.11 11.99 5.76 4.11 
1000 5.28 4.91 4.88 11.41 8.98 5:21 9.35 20.89 13.87 6.36 4.46 
2000 5.66 5.28 DS 1272 9.88 527) 10.3 23.68 15.75 6.96 4.81 
5000 6.16 5.77 5.93 14.41 11.06 6.35 11.5 27.33 18.23 7.75 5.26 
10 000 6.54 6.16 6.37 UDI 1122 6.83 12.41 30.22 20.07 8.38 5.63 
ro RH 26 Pánuco RH 26 Valle de México RH 27 
G1 G2 G3 Gl G2 G3 G2 G3 G1 G2 G3 
Años 
2 0.84 0.79 0.81 0.91 0.81 0.8 0.88 0.86 0.83 0.85 0.77 
5 ME 1.47 1.38 1.42 1.41 1.71 Al 1.6 1.28 1.4 127 
10 1.9 2.05 2.01 1.76 2.08 2.22 1.93 2.08 1.78 1.88 1.84 
20 2.42 2.79 2.67 2.08 2.6 2.66 2.32 2.54 2.48 2.23 2.91 
50 2.96 3.99 3.29 2.5 3.13 3.19 2.84 3.15 3.2 2.61 3.92 
100 3.33 4.92 8%) 2.82 3.49 3.58 3.23 3.6 3.67 2.88 4.57 
200 3.68 5.84 4.09 3.13 3.84 3.96 3.61 4.04 4.12 3.14 5.18 
500 4.14 7.02 4.6 3.54 4.3 4.47 4.12 4.64 4.71 3.49 5.97, 
1000 4.49 7.92 4.98 3.86 4.64 4.85 4.5 5.08 5.14 3.74 6.56 
2 000 4.84 8.8 D0:3Z 4.17 4.98 5.23 4.89 5.53 5.58 4 7.15 5 
5000 5.29 9.99 5.86 4.58 5.45 5.73 5.4 6.12 6.15 4.35 7.93 2 
10000 5.63 10.86 6.24 4.9 5.79 6.14 5.78 6.57 6.58 4.59 8.54 E 
S 
e RH 28 RH 30 RH 36 y 37 a 
G1 G2 G3 G1 G2 G3 G1 G2 G2 E 
Años Q (m/s) | Q (m/s) | O (m/s) | O (m/s) | Q (m/s) | Q (m/s) | Q (m/s) | Q (m/s) Q (m/s) o 
2 0.89 0.92 0.76 0.9 0.87 0.77 0.72 1.02 1.01 E 
5 eS 1.37 1.17 1.21 1127 11,222 1.61 1.75 1.8 E 
10 1.72 1.67 1.45 1.49 1.61 1.6 2.32 2.24 2.32 be 
20 1.96 1.96 JA 1.84 2.07 2.32 2.88 2 JA 2.81 E 
50 2.19 2.33 2.05 2.34 2.85 5.66 3.54 3.32 3.46 e 
100 2.35 2.61 2.31 27 3.44 8.18 4.02 3.77 3.94 > 
200 2.51 2.89 2.56 3.04 4.02 10.57 45 4.23 4.42 = 
500 DTO! 3.26 2.9 3.49 4.77 13.63 5.12 4.83 5.06 El 
1 000 2.89 3.53 3.15 3.83 5.33 15.93 5.59 5.28 5.54 = 
2 000 3.05 3.81 3.41 4.17 5.89 18.19 6.05 US 6.02 3 
5000 3.28 4.18 3.74 4.61 6.63 21.28 6.68 6.33 6.65 3 
10 000 3.45 4.45 3.99 4.94 UE2a 23.38 Ava 6.78 7.13 S 
É 
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delos regionales para estimar la MOMIA en 
función de dichas características (cuadro 5). En 
el cuadro 6 se presenta el coeficiente de deter- 
minación obtenido con cada modelo. 

En los cuadros 5 y 6, la ecuación potencial (3) 
se logró usando una hoja de cálculo ExcelO al 
dibujar la MQMIA de las hidrométricas de cada 
grupo de regiones contra el área de la cuenca 
correspondiente, y proponiendo una línea de 
tendencia potencial de la que se estableció su 
ecuación y coeficiente de determinación. El 
siguiente modelo fue del tipo de la ecuación 
(4), cuyos parámetros se definieron usando 
algoritmos genéticos; dicho modelo, además 
del área de la cuenca, considera el tiempo de 
concentración en la precisión de la MOMIA. 


En las ecuaciones del tipo (5) se consideró el 
volumen medio llovido (al tomar en cuenta el 
área de la cuenca y precipitación promedio en 
la cuenca, calculada con el promedio espacial y 
la media de las precipitaciones diarias máximas 
anuales, además del tiempo de concentración y 
efectos de infiltración en la cuenca al añadir al 
espesor de suelo S). 

Cabe destacar que en el tercer modelo se 
agrupó el área con la precipitación y con ello se 
logró mejorar los coeficientes de determinación, 
aumentando sólo en 1 el número de parámetros 
que se requiere calcular. 

En el caso de cuencas no aforadas y depen- 
diendo de las características fisiográficas y los 
datos de precipitación, primero se debe estimar 


Cuadro 5. Ecuaciones para estimar la MOMIA usando un modelo potencial y modelos no lineales a partir de AG. 


RH Ecuación potencial 


1y3 | MOMIA= 0.0002*A17%5 


Ecuación AG con A y tc 


MOMIA az 0.30529* 4079293 ¿-0.09583 


Ecuación AG con V, tc y S 


MOQMIA E 0.0001 Y 0-31044.0.7405G3.0198 


7y8 MQMIA = 142.39* 40005 


MOMIA = 0.30529* 41:60802+ ¿ ¿-2.51320 


MOMIA = 877.7123V1-1838fg-1.3504 G-4.4633 


9 MOMIA = 11.66* 400 


MOQMIA pan 2.44158* A1-0198% p -1.03305 


MOMIA = 1 .2208 Y 12204f -0.515-3.5978 


10 MOMIA = 20.745* 40412 


MOMIA = 6.71417* A073861%(-052589 


MOMIA = 0.6105 V07187f(-0.35865-0.4141 


11 MOMIA = 5.6048* 403886 


MOQMIA = 0.61047* A1-08572* -2.10456 


MOMIA = 0.3053 V 1-8204f-1.3797G-4.6715 


12 MOMIA = 0.1459*408%6 


MOMIA =0.30529* A!51708%p(-173559 


MOMIA = 0.3053 Y 15354 1.3102 G-3.5441 


13 a16 | MQMIA =1.3419*407577 


MOMIA a 8.54528* 4059389 ¿¿-0.18403 


MOMIA a 0.3053V153354.fc-0-1883G0.1273 


19 MOMIA = 2.0168* 408158 


MOMIA = 0.30529*A123254*p-036105 


MOMIA = 0.0001 Y 1-256440-187160.4636 


18 MOMIA = 26.475* 407% 


20 y 21 | MOMIA = 11.974* 404763 


MOMIA = 0.61047* A12694%(-127964 


MOMIA = 36.01194* 4020475052705 


MOMIA a 3.3571 Y 0.16304-0-1022G-0.1999 


MOMIA = 0.001 Y 0-91654 0.3745 G2.551 


22 MOMIA = 27.621* 4096 


MOMIA pa 38.45342* 4044545* (-0.52248 


MOMIA E 0. 61 05 y 0.5587 0.5667 C0.4416 


28) MOMIA = 6.442* 406202 


MOMIA = 3.35714* A0829*p-036105 


MOMIA = 46.083 V/0-0488(-0.457560.1926 


24 MOMIA = 0.256* 40801 


MOMIA = 0.30531* 415973 (-15288 


MOMIA = 0.3053 V1-2622f(-1.47255-1.7765 


20 MOMIA = 115.98* 413% 


MOMIA = 107.7304*401473*p(0 


MOMIA = 9 974.3645 V%tco24115-1.7038 


26P MOMIA = 33.527* 403732 


MOMIA = 7.3245* 409372 f(-1.0324 


MOMIA == 10.9868 Y 0.8591 f-0.7047 C-1.8973 


26 VM | MQMIA =5.1785*407% 


MOMIA = 11.2919*A00162*p(1.0361 


MOMIA = 7.0194 V%fc09488 502963 


27 MOMIA = 1.6791* 476% 


MOMIA = 31.7394*A02689%(0.4324 


MOMIA = 0.3053 V0445fc0-414150.5802 


28 y 29 | MOMIA = 69.115* 403% 


MOMIA = 60.73193* A09651%p(0 


MOMIA = 0.6105 Y 0.492 0.276860.1541 


30 MOMIA = 40.932* 409% 


MOMIA = 0.61047* A081301%¿(0.13215 


MOMIA a 3.6623 Y 0-8805 (0.1621 5-2.9502 


36 y 37 | MOMIA = 18.941*403%20 


MOMIA = 1.83121* 4093387*p072725 


MOMIA = 10 000 V03316f-0-1279 5-3.6222 


Nota: VM: Valle de México; P: Pánuco. 
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Cuadro 6. Coeficiente de determinación de los modelos del cuadro 5. 


RH Ecuación potencial Ecuación AG con A y tc Ecuación AG con V, tc y S 
1y3 0.442 0.544 0.897 
7y8 0.038 0.997 1 
9 0.684 0.749 0.994 
10 0.661 0.697 0.690 
11 0.172 0.769 0.903 
12 0.769 0.928 0.976 
13a16 0.663 0.703 0.810 
19 0.730 0.760 0.822 
18 0.411 0.817 0.832 
20 y 21 0.390 0.443 0.827 
22 0.115 0.402 0.524 
23 0.488 0.504 0.513 
24 0.436 0.648 0.919 
25 0.522 0.528 0.629 
26 P 0.630 0.736 0.870 
26 VM 0.364 0.647 0.652 
27 0.259 0.365 0.454 
28 y 29 0.318 0.325 0.438 
30 0.571 0.864 0.940 
36 y 37 0.430 0.790 0.984 


Nota: VM: Valle de México; P: Pánuco. 


la MOMIA utilizando las ecuaciones indicadas 
en el cuadro 5, y después se utilizaría el factor 
correspondiente a la región y grupo del sitio 
analizado. 

En el cuadro 6 se observa que el modelo de 
la ecuación potencial, que sólo toma en cuenta 
el área, da coeficientes de determinación bajos, 
menores a 0.2, en las regiones agrupadas 7 y 8, y 
en las regiones hidrológicas 11 y 22; esto es, que 
en esas regiones, los parámetros tc, V y S deben 
considerarse en la estimación de la MOMIA. 


Ejemplo de aplicación de la ecuación regional 
para la MQMIA en la RH19 


En 2013 se registraron eventos de precipitación 
extraordinarios en Guerrero, en particular 
en el río La Sabana, que ocasionaron fuertes 
caudales. Si se toma en cuenta la ecuación re- 
gional obtenida en este estudio, la estimación 
de la MOMIA para la estación 19018 Tuncingo 


tendría una magnitud aproximada de 325 m?*/s 
(cuadro 7). Dicho valor se calcula al considerar 
que dicha hidrométrica se encuentre en la región 
hidrológica 19, cuya ecuación regional (figura 
7) estima la MOMIA, en m/s, en función del 
volumen llovido en km?*.mm, obtenido como el 
producto del área de la cuenca drenada hasta 
el sitio de aforo, en km”, por la media de las 
precipitaciones diarias máximas anuales hp 
en mm, el tiempo de concentración en horas y 
espesor del suelo S en cm. 

Al multiplicar la MOMIA de 325 m?/s por el 
valor del gasto máximo instantáneo normaliza- 
do de diseño para periodos de retorno de 50 y 
100 años, obtenido a partir del análisis regional 
(ver cuadro 4), se determinan gastos máximos 
instantáneos de 1 335.70 m*/s y 1 621.69 m?*/s, 
para los periodos de retorno de 50 y 100 años, 
respectivamente (cuadro 8). 

Por otra parte, la Conagua (2013), utilizando 
el registro histórico de gastos máximos anuales 
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Medida calculada vs. media 19 


1200 


MOMIA medida, m3/s 
Ha [e] 00 a 
2.8 8 8 


155] 
e 
[a] 


0 200 400 600 


800 1000 


y =0.937x + 42.457 
R2=0.8219 


+ Media calculada 
US. 
media medida 


— Identidad 


— Lineal (media calculada 
1200 vs. media medida) 


Media de Qmi calculada, m3/s 


Figura 7. Ecuación regional para la RH19. 


Cuadro 7. MOMIA obtenidos con la ecuación regional para la estación 19018 Tuncingo, Guerrero. 


MOMIA con 
Área hp VA tc Ss ec. regional | Grupo RH13 
1 
avs Nombre km? hp mm km?.Mm h espesor AG al6 y 19 
m/s 
19018 Tuncingo 399 127.49 50 868.51 5.44 10 324.99 3 
Cuadro 8. Gastos máximos instantáneos obtenidos con análisis regional para la estación 19018 Tuncingo. 
Ajustes de MOMIA normalizados O de diseño m?*/s 
Clave Nombre Tr Tr 
50 años 100 años 50 años 100 años 
19018 Tuncingo 4.11 4.99 1 335.7 1 621.69 


de 1969 a 2013 (este último dato estimado con 
la huella del agua registrado por la creciente 
de septiembre de 2013 ante los eventos de los 
huracanes Manuel e Ingrid), determinó gastos de 
diseño de 1 207 m*/s para el periodo de retorno 
de 50 años y de 1 389 m?/s para el periodo de re- 
torno de 100 años; es decir, valores ligeramente 
menores a los calculados en el cuadro 8 con el 
estudio regional. 

Además, se realizó una comparación utili- 
zando un modelo de pronóstico de escurrimien- 
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to (MPE) (Domínguez et al., 2008) para estimar 
el gasto máximo en el río La Sabana, a la altura 
de la estación hidrométrica 19018, con el evento 
de precipitación total de 619.33 mm ocurrido 
entre el 12 y 18 de septiembre de 2013 (cuadro 
9 y figura 8). Al considerar los parámetros del 
cuadro 10, se obtuvo el hidrograma que se 
muestra en la figura 9, con un gasto pico de 
1 419.42 m*/s, el cual corresponde al gasto máxi- 
mo instantáneo comprendido entre los periodos 
de retorno de 50 y 100 años. 
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Cuadro 9. Precipitación y gastos obtenidos en el MPE. 


19018 Tuncingo, cuenca del río La Sabana 


Fecha hp (mm) Q total(m//s) 
12/09/2013 14:00 0 0 
12/09/2013 20:00 17.26 0 
13/09/2013 02:00 2.63 Si02 
13/09/2013 08:00 0.12 5.19 
13/09/2013 14:00 4.72 1.91 
13/09/2013 20:00 1.16 2.88 
14/09/2013 02:00 32.5 3.34 
14/09/2013 08:00 35.6 61.82 
14/09/2013 14:00 35.75 200.43 
14/09/2013 20:00 81.1 331.45 
15/09/2013 02:00 82109) 714.29 
15/09/2013 08:00 87.72 776.19 
15/09/2013 14:00 101.15 882.25 
15/09/2013 20:00 49.7 1 419.42 
16/09/2013 02:00 40.69 1 202.35 
16/09/2013 08:00 20.18 688.91 
16/09/2013 14:00 24.62 366.67 
16/09/2013 20:00 38.55 190.39 
17/09/2013 02:00 3.59 312.03 
17/09/2013 08:00 9.5 240.62 
17/09/2013 14:00 0.2 34.37 
17/09/2013 20:00 4.91 
18/09/2013 02:00 0.7 


Cuadro 10. Parámetros de calibración. 


Parámetros de calibración Río La Sabana 
AT (h) 6 
% 0.2 
Fs 1 
Fo 0.97 
tc (h) 5 


K (h) 


Conclusiones y recomendaciones 


En este trabajo se utilizó la media de los gastos 
máximos anuales (MOMIA) como concepto 
eje, de tal manera que al relacionarla con 
características medibles en cualquier cuenca 


(área, tiempo de concentración, promedio 
espacial de la media de precipitaciones máxi- 
mas anuales, retención potencial máxima) es 
posible estimarla de forma sencilla, aun para 
cuencas no aforadas y, por otro lado, al multi- 
plicarla por los factores indicados en el cuadro 
4, se obtiene una estimación de los máximos 
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Precipitación en la cuenca 
La Sabana 


Figura 8. Precipitación en la cuenca La Sabana. 


Hidrograma directo 
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Figura 9. Hidrograma obtenido a partir de la precipitación. 


instantáneos asociados con distintos periodos 
de retorno. 

Los resultados relativos a la estimación de 
la MOMIA se resumen en los cuadros 5 y 6. 
En el primero se presentan las ecuaciones que 
permiten estimarla a partir de características fá- 
cilmente medibles en cuencas no aforadas, con- 
siderando tres posibles modelos de regresión. 


ISSN 0187-8336 


La confiabilidad de los modelos se refleja en 
el coeficiente de determinación, aunque también 
es importante buscar su simplicidad y lógica de 
las variables involucradas, de tal manera que 
un análisis del cuadro 6 muestra que el modelo 
más simple (en función sólo del área de cuenca) 
podría aplicarse en los casos de las regiones 9, 
10, 12 y 19; por otra parte, el tercer modelo (con 
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V, tc y S), que involucra precipitación, es el me- 
jor en todos los casos, excepto para las regiones 
7, 8, 10, 18, 23 y 26, en las que el incremento en 
el coeficiente de determinación es pequeño. 

Por otra parte, para varios grupos de 
regiones se obtuvieron exponentes positivos 
asociados con la retención potencial máxima 
y/o tiempo de concentración, lo cual no parece 
lógico, por lo que se recomienda estudiar posi- 
bilidades, como la de acotar los valores de estos 
exponentes al usar la optimación vía algoritmos 
genéticos, así como la influencia en resultados 
de la correlación entre el área de las cuencas y 
su tiempo de concentración. 

Para la estimación de los coeficientes regio- 
nales mostrados en el cuadro 4 se utilizaron 305 
estaciones asociadas con cuencas “vírgenes” 
(es decir, no afectadas de modo sensible por 
zonas de regulación naturales o artificiales) y 
que contaran con más de 20 años de registro; 
sin embargo, los resultados pueden usarse para 
estaciones con menos años de registro (p. ej., del 
orden de 15 años), debido a que en esos casos 
no sería necesario hacer un análisis estadístico 
que requiera momentos de orden mayor a uno, 
y sólo se necesitaría estimar la media de los 
valores máximos anuales registrados. 

Se recomienda actualizar el estudio de for- 
ma periódica (p. ej., cada cinco años, como se 
sugiere para la revisión de avenidas de diseño). 

Asimismo, con los datos generados en este 
estudio (se cuenta con bases de datos para las 
más de 300 cuencas analizadas, que contienen 
valores depurados de los gastos máximos 
anuales instantáneos y para duraciones de l a 
60 días, así como sus áreas, también depuradas, 
tiempos de concentración, numero de curva, 
etc.) pueden realizarse otros (p. ej., los autores 
de este trabajo están llevando a cabo un estudio 
estadístico de los coeficientes de escurrimien- 
to). Para los investigadores interesados en el 
tema, se cuenta también con valores de lluvia 
máximos anuales para duraciones de 1 a 60 días 
depurados para más de dos mil estaciones, y 
las extrapolaciones correspondientes a diversos 
periodos de retorno y duración. 
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